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We present a study of recent solar neutrino data using 
a Bayesian method. Assuming that only v e are observed in 
the Super-Kamiokande experiment our results show a marked 
supression of the survival probability at about 1 MeV, in good 
agreement with x 2_ based analyses. When the detection of 
f M by Super-Kamiokande is taken into account, assuming v e 
to i/ M oscillations, we find the largest suppression in survival 
probability at about 8.5 MeV. 

PACS numbers : 26.65.+t,13.15.+g 



One of the most intriguing problems of the past two 
decades has been the observation of a deficit of neutrinos 
of solar origin as compared to the predictions of stan- 
dard solar models Jl]-|j| . Many attempts have been made 
to explain this discrepancy either as a consequence of 
astrophysical processes or new physics such as neutrino 
oscillations in vacuum Q] or in matter |^| . The astrophys- 
ical solutions have not been successful ||. However, the 
solutions that invoke new physics provide excellent de- 
scriptions of the solar neutrino data In this paper 
we answer the following question: how well do we know 
the neutrino survival probability as a function of the neu- 
trino energy? The answer is relevant because it tells us 
in what way, and to what degree, the solar neutrino data 
constrain the models that seek to explain the neutrino 
deficit. We answer the question using a Bayesian method. 
Our perspective here is broader than that of Ref. || in 
which a Bayesian method was used to analyze the MSW 
model. 

In our analysis we assume that the solar neutrino 
spectrum is that predicted by the standard solar mod- 
els ||i|fl[L0[ . However, it is known that the spectrum is 
insensitive to the details of these models JllJ] . The ex- 
perimental data are from Homestake (CI) fll2|, SAGE 
(Ga) |l^_GALLEX (Ga) Q and Super-Kamiokande 
(H 2 0) 



15|jlq]. These results together with the predic- 
tions of the standard solar model of Bahcall and Pinson- 
neault jlj are shown in Table I. Our method of analysis 
does not require the imposition of the solar luminosity 
constraint. In accordance with our minimalist approach 
we choose not to impose it. 

Bayes' theorem, P(H\D,I) = C{D\H,I) P{H\I) 
j J C(D\H,I) P(H\I), gives a prescription for calculat- 
ing the posterior probability P(H\D, I) of an hypothesis 
H, given measured quantities D and prior information /; 



L is the likelihood function assigned to D and P(H\I) is 
the prior probability assigned to H. The integration in 
the denominator is over all hypotheses of interest. 

The solar neutrino rate Si for the Chlorine and Gallium 
experiments is given by, 



J E, 



(1) 



where $j is the total flux from neutrino source j, 4>j is 
the corresponding normalized neutrino energy spectrum, 
cr, is the cross-section for the ith experiment, E t hi is its 
threshold energy (see Table I) and V{E U ) is the neutrino 
survival probability. 

For the Super-Kamiokande experiment we use their re- 
ported measurement of the electron recoil spectrum pro- 
duced by the 8 P> neutrinos that spans the range 
6.5 to 20 MeV. If the v e deficit is caused by v e oscil- 
lations to tV T ) then one must take account of the fact 
that Super-Kamiokande is sensitive to (but does not dis- 
tinguish between) all flavors of neutrino. On the other 
hand, if the v e disappear through a mechanism that does 
not result in other detectable particles, for example by 
oscillating into sterile neutrinos, then the measured rate 
is to be ascribed to the v e flux only. We consider both 
possibilities. 

The measured electron recoil spectrum N(T) is given 

by 



N(T) = N 



dT'R(T\T') 



(2) 



E" 



< I dE v <t> B (E u )[V(E v )a e (T',E v ) 
+ (l-V(E I/ ))a fl (T',E u )}, 

where R(T\T') is the Super-Kamiokande resolution func- 
tion (which can be approximated by a Gaussian with 
mean V and standard deviation LSy^T'/lOMeV) Uj), 
T = E e — m e and T' are the measured and true elec- 
tron kinetic energies, respectively, with E e and m e the 
electron energy and mass. The quantity 4>b is the nor- 
malized neutrino energy spectrum from the S P> reaction, 
and <r e and are the v e and Vn differential electron 
scattering cross-sections respectively. Given a neu- 
trino energy E v the electron can assume a maximum 
kinetic energy of T max '(E u ) = 2E*/(2E V + m e ), while 
the minimum neutrino energy for a fixed T' is given by 



1 



TABLE I. Measured and predicted solar neutrino capture 
rates [16], in units of SNU (=10 -36 captures/atom s _1 ), for 
the radiochemical experiments, and neutrino fluxes, in units 
of 10 6 cm -2 s _1 , for Super-Kamiokande. 



Experiment 


E t h 


Flux Rates 




(MeV) 


Measured 


Predicted 


Homestake 


0.814 


2.54±0.16±0.14 


9.30±1.30 


K+ 37 C1^ e-+ 37 Ar) 








GALLEX 


0.233 


76.2±6.5±5 


137±8 


SAGE 




73±8.5±H 


If 


(i/ e + 71 Ga -> e~+ 71 Ge) 






Super-Kamiokande 


6.5 


2.44±0.06+o;o9 




(v e +e~ — > e~+v e ) 









E™ in (T') = [T' + y/{T'(T' + 2m e ))]/2] E™ ax is the max- 
imum neutrino energy, which we take to be 20 MeV. The 
constant iVo is a normalization factor that depends on 
which units are used for the event rate. In writing Eq. (2) 
we have assumed v e to oscillations. If the measured 
flux is due to v e only the recoil spectrum is given by 
Eq. (2) with the term proportional to omitted. The 
event rate Si in the ith electron energy bin is simply 
N(T) integrated over that bin. 

We note that each experiment is sensitive to differ- 
ent parts of the neutrino energy spectrum. This is ev- 
ident from the (normalized) plots of cr(E v ]Y^^j4>j(E 1 y) 

(or ME,) x CHI dT Io iEv) dT'R(T\T') a e (T',E u ) 
for Super-Kamiokande) shown in Fig. [|. We also note 
the existence of regions where the spectral sensitivity is 
essentially zero. 

The likelihood function C(D\H,I) is assumed to be 
of Gaussian form g(D\S, £), where D = (Di, . . . , Dig) 
represents the 18 data — 2 rates from the radiochemi- 
cal experiments plus 16 rates from the binned Super- 
Kamiokande electron recoil spectrum, £ is the 18 x 
18 error matrix for the experimental data and S = 
(Si, . . . , Sxs) represents the predicted rates. The error 
matrix is deduced from the data given in Ref. JT6|. We 
take the prior probability to be constant. 

We extract the survival probability using two different 
methods: binned and parametric. In the binned method 
we divide the neutrino energy spectrum into 12 bins be- 
tween 0.2 and 20 MeV, chosen so that the survival prob- 
ability within each bin is approximately constant. With 
a minor algebraic manipulation of Eq. (1), the neutrino 
rates Si and 52 for the Chlorine and Gallium experi- 
ments, respectively, reduce to a weighted sum of Vk (the 
survival probability in bin k): 



Si ~J2 Vk Y. $ J / ViWh {E v )dE v , (3) 



while for Super-Kamiokande, assuming v e to oscilla- 
tions, the neutrino rate is expressed as sum of the sixteen 
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1. Spectral sensitivity as a function of the neutrino 



spectral values S3 to Sis, each containing two terms: one 
with (T^ and the second with the difference between a e 
and (7 M (see Eq. (2)). If electron neutrinos are the only 
detectable particles then one gets a similar expression 
with the <7 M terms omitted. 

From Bayes' theorem we compute the posterior prob- 
ability P(V\D, I), that is, the probability of the hypoth- 
esis H that the set of parameters V = (Vi, ■ ■ ■ , V12) 
have specified values. For the Gallium rate we use the 
weighted average of the SAGE and GALLEX results. 
The posterior probability P(V k \D,I) for each Vk is ob- 
tained by marginalizing (that is, integrating) P(V\D,I) 
over the remaining Vk- We use the mean of P(Vk\D, I) 
as our best estimate of Vk- The results, assuming that 
the measured rate is produced by v e only, are shown 
in Fig. ||(a). We see that in the first two bins, that 
span the pp neutrino spectrum, the survival probability 
is close to unity. There is a marked suppression in the 
bin (0.8-1.5 MeV) that contains the main 7 Be line and 
a moderate suppression in the 8 B> spectrum (bins span- 
ning 6-12 MeV). These results agree with previous analy- 
ses |?],|8| . But our result provides additional information, 
namely: the precision with which the survival probability 
is known — as a function of the neutrino energy, irrespec- 
tive of the precise origin of the neutrino deficit. Moreover, 
since our binned method makes no assumption about the 
form of the survival probability it makes no unwarranted 
inferences. Where there is little spectral sensitivity (in 
the intervals 0.4-0.8 MeV and 1.5-4.5 MeV) our method 
infers correctly that little information can be extracted. 

If, however, we assume oscillations into active neu- 
trinos (vy) the inferred form of the survival probability 
changes. Figure 2(b) shows the extracted survival proba- 
bility for such a case. Here, we see a marked suppression 
at a higher neutrino energy (i.e., at around 8.5 MeV). 

In the parametric method, which is motivated by the 
conclusions of previous analyses |7P,ll8|], we write the 
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FIG. 2. Our best estimate of the survival probability (dark 
dashed line) vs neutrino energy for the binned method assum- 
ing (a) the neutrino flux consists of v e only and (b) i/ e plus 
Vfi, due to oscillations. The 68.3% (shaded area) and 90% 
(solid line) confidence intervals, about our best estimate, are 
also shown. 

survival probability as a sum of two finite Fourier series: 

7 

V(E„\a) =Y, a r+icos(r%E v /Lx) (4) 
/ (l + exppk-£i)/6]) 

3 

+ a r+ QCOs{ri:E v / L 2 ). 

r=Q 

The first term in Eq. (4) is defined in the interval 0.0 
to L\ MeV — and suppressed beyond L\ by the exponen- 
tial, while the second term covers the interval 0.0 to L 2 
MeV. We divide the function in this way so that it can 
accomodate a survival probability that varies rapidly in 
the interval 0.0 to L%. Holding the parameters L%, L 2 
and b fixed at 1.0, 15.0 and 0.1 MeV, respectively, we 
compute the posterior probability P{a\D 1 I) for the pa- 
rameters a = (oi, . . . , 012). The data are the same as 
those used for the binned method; but instead of a linear 
sum in Vk we have a linear sum in the parameters a. 

The theoretical uncertainties, which so far we have 
neglected, can be incorporated by treating the fluxes 
$ = ($1,...) as parameters with an associated prior 
probability, P($\I), that encodes the flux predictions 
and whatever correlations exist amongst them. We rep- 
resent our knowledge of the fluxes by a multivariate 
Gaussian prior probability P($\I) = X^), where 

$° = ($5, ■ ■ ■ j ®&) is the vector of flux predictions and 
is the corresponding error matrix [^9| . For simplicity 
we neglect correlations; hence is diagonal. 

The posterior probability (using a constant prior 
probability for the parameters a is now given by 
P(a,$\D,I) = C(D\a,$,I) P(*|J) / J* 0) $ C(D\a, I) 
P($>\I), which when marginalized with respect to $ gives 




1 10 

Neutrino Energy (MeV) 



FIG. 3. Survival probability vs neutrino energy for the 
parametric method assuming (a) the neutrino flux consists of 
v e only and (b) v e to v^. All uncertainties are included. See 
Fig. ^ for other details. 

P(a\D,I). The survival probability is estimated using 
P(E U \D,I) = $ a V{E y \a) P(a\D,I). Figure |(a) shows 
the survival probability and its uncertainty, which now 
includes both the experimental and theoretical errors. 
Again, the general form of the survival probability ob- 
tained, assuming the neutrinos detected consist of v e 
only, agrees with the inferences from previous analy- 
ses and our binned method. But again, unlike previ- 
ous analyses we have detailed information about the sur- 
vival probability and its uncertainties as a function of 
the neutrino energy. Figure [|(b) shows our results, for 
the survival probability, assuming oscillations into active 
neutrinos. 

One useful check of the reliability of both calculations, 
and the flexibility of the function in Eq. (4) , is to set the 
measurements D equal to their predicted values and to 
verify that the extracted survival probability is unity in 
the energy ranges where data exist and is 0.5 with an 
error of ±0.3 where there are no data. Both calculations 
have been verified successfully in this way. We note the 
general agreement between the binned and parametric 
methods in the energy regions with adequate spectral 
sensitivity. 

The calculations presented here are based on data from 
1997. Our preliminary calculations using the data pre- 
sented at the recent Neutrino '98 conference |20| , and 
the 1998 Bahcall and Pinsonneault predictions |2l| , yield 
similar conclusions, as illustrated in Fig. ||. The Bayesian 
method offers a well-founded way to compute the relative 
probabilities of different models of the survival probabil- 
ity. We plan to provide such calculations in a longer 
paper. 

In summary, we have used recent solar neutrino data 
and standard solar model predictions to extract the neu- 
trino survival probability and its uncertainty as a func- 
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FIG. 4. Survival probability vs neutrino energy for the 
parametric method using 1998 data and solar model[21]. All 
uncertainties are included. See Fig. 3 for other details. 

tion of the neutrino energy under two broad alternate 
assumptions: 1) the measured flux is comprised of v e 
only or 2) it is a mixture of active neutrinos arising from 
v e oscillations. Under assumption 1), we find that the 
survival probability is most precisely determined around 
0.3 MeV and has a value of 0.79±0.13 (and 0.28±0.03 at 
8.3 MeV) with experimental and theoretical uncertain- 
ties included. However, the uncertainties elsewhere are 
considerably larger. Our results suggest that the data 
could accomodate a wider range of models than hitherto 
have been considered. In particular, the fact that the 
MSW model provides a good description of the data is 
clearly a strong point in its favor; but it is at present not 
a decisive one. We therefore eagerly await further results 
from Super-Kamiokande and the first results from SNO 
II- 
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